DUDLEY WOX LIBRARY 
NAVAL POSTGRADUATE SCHOOL 
MONTEREY, CA 93940 



NAVAL 



POSTGRADUATE SCHOOL 

Monterey, California 




THESIS 



SELF GENERATED MAGNETIC FIELDS 
IN LASER- PRODUCED SHOCK WAVES 

by 



Cem Yrldmm Parlar 
December 1978 



Thesis Advisor: 



F . Schwirzke 



Approved for public release; distribution unlimited. 



T 18739 



UNCLASS TFT Fn 



SECURITY CLASSIFICATION OF THIS PACE (Whon Dmto Entorod) 



REPORT DOCUMENTATION PAGE 


READ INSTRUCTIONS 
BEFORE COMPLETING FORM 


1. REPORT NUMBER 


2, GOVT ACCESSION NO. 


3. RECIPIENT'S CATALOG NUMBER 


4. TITLE (ltd Submit) 

Self Generated Magnetic Fields in 
Laser-Produced Shock Waves 


3. TYPE OF REPORT A PERlOO COVERED 

Master's Thesis; 
December 1978 


«. PERFORMING ORG. REPORT NUMBER 


7. AUTNORf»> 

Cem Yildinm Parlar 


8. contract OR GRANT NUMBER^*) 


9. PERFORMING ORGANIZATION NAME ANO AOORESS 

Naval Postgraduate School 
Monterey, California 93940 


10. PROGRAM ELEMENT, PROJECT, TASK 
AREA 6 WORK UNIT NUMBERS 


1 1. CONTROLLING OFFICE NAME ANO AOORESS 

Naval Postgraduate School 
Monterey, California 93940 


12. REPORT DATE 

December 1978 


<3. NUMBER OF PAGES 

107 


14. MONITORING AGENCY name a ADDRESS*’// diiimront from Controlling Oftico) 


IS. SECURITY CLASS. (o( this rtpanj 

Unclassified 


!5«. OECLASSI F| CATION/ DOWNGRADING 
SCHEDULE 


16. DISTRIBUTION STATEMENT (of thio Ropoet) 

Approved for public release; distribution unlimited. 

c 

• 


17. DISTRIBUTION STATEMENT { of tho obotroct muorod tn Block 20, // dlfforont from Report) 


18. SUPPLEMENTARY NOTES 


18. KEY WORDS (Continue on towotoo sido it nm coooory and i don t tty by block numb or) 

Laser-Produced Shock Waves 
Self-Generated Magnetic Fields 


20. ABSTRACT (Continue on rovoewm wide it nocooooey and i don t tty by block number) 

The generation of spontaneous magnetic fields in a plasma 
produced by a 5.0 Joule, 22 nanosecond laser pulse and the 
generation of the reversed magnetic fields after several hundred 
nanoseconds due to the interaction between the laser produced 
plasma and a background plasma was investigated. 

The generation of the reversed magnetic fields coincides with 
the steepening of the expanding plasma front. This steeoenina 



DD I jan *73 1473 EDITION or I MOV •• !S OBSOLETE 
S/N 010J-014-4S01 | 



UNCLASSIFIED 



1 



SECURITY CLASSIFICATION OF THIS FACE (Whit Data Mnfrii) 



UNCLASSIFIED 



{aeuwTv cl A ttiFic* tion or tmis wtotrw*— > n»«« 



ABSTRACT (Cont'd) 

is caused by a collisional shock wave which is formed during 
the interaction of the expanding laser-produced plasma with 
the photoionized background gas. The analysis of the source 
term for magnetic field production indicates that the genera- 
tion of the reverse fields is due to an axial temperature 
gradient and a radial density gradient. This axial tempera- 
ture gradient is a result of the heating of the plasma in 
the collisional shock wave. 

The width of the shock was measured to be 0.2 cm. The 
magnitude of the maximum reverse field was measured at 5.0 
Torr (H ) background pressure to be 150 Gauss. 



DD Form 1473 




UNCLASSIFIED 



2 *teu»M*V CUAMIUCATION 0^ T*l» Pratt**** £«<•'•«*) 



Approved for public release; distribution unlimited. 

Self Generated Magnetic Fields 
in Laser-Produced Shock Waves 

by 

Cem Yildirim Parlar 
Lieutenant, Turkish Navy 
B.S., Naval Postgraduate School, 1978 



Submitted in partial fulfillment of the 
requirements for the degree of 



MASTER OF SCIENCE IN PHYSICS 



from the 



NAVAL POSTGRADUATE SCHOO 
December 1978 



9 \%^ 

' - \ 








DUDLEY KNOX LIBRARY 

naval postgraduate school 

MONTEREY. CA 93940 



ABSTRACT 



The generation of spontaneous magnetic fields in a 
plasma produced by a 5.0 Joule, 22 nanosecond laser pulse 
and the generation of the reversed magnetic fields after 
several hundred nanoseconds due to the interaction between 
the laser produced plasma and a background plasma was 
investigated . 

The generation of the reversed magnetic fields coincides 
with the steepening of the expanding plasma front. This 
steepening is caused by a collisional shock wave which is 
formed during the interaction of the expanding laser- 
produced plasma with the photoionized background gas. The 
analysis of the source term for magnetic field production 
indicates that the generation of the reverse fields is due 
to an axial temperature gradient and a radial density 
gradient. This axial temperature gradient is a result of 
the heating of the plasma in the collisional shock wave. 

The width of the shock was measured to be 0.2 cm. The 
magnitude of the maximum reverse field was measured at 5.0 
Torr (H^) background pressure to be 150 Gauss. 
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I. 



INTRODUCTION 



A dense, energetic plasma can be produced by focusing 
the pulse from a powerful laser onto a solid target. Asso- 
ciated with this laser-produced plasma (hereafter referred 
to as laser plasma) spontaneously generated magnetic fields 
have been reported and investigated [Refs. 22,24]. These 
investigations imply that spontaneously generated magnetic 
fields arise generally during the production of a laser 
plasma. The strength of these fields may vary with target 
geometry [Ref. 31], background gas pressure [Ref. 3] and 
amount of energy delivered to the target by the laser. In 
a computer simulation experiment, Widner [Ref. 30] has 
observed a spontaneous field intensity of the order of 1 
Megagauss . 

Stamper et. al. [Ref. 25] have reported spontaneous 
magnetic fields, of the order of 0.6 Megagauss, generated 
in the plasma near the focus of a high-irradiance Nd-laser 
beam. They used a Raman-shifted, three channel Faraday- 
rotation diagnostic system and observed that the spontan- 
eously produced magnetic fields depend on a variety of 
experimental parameters, such as target material and 
geometry, laser irradiance and polarization, focal position, 
and prepulse effect. 

The influence of magnetic fields on laser plasmas is an 
important topic in laser-fusion research. The presence of 
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such fields in a plasma has effects on electronic heat 
conduction and the generation of plasma instabilities. 

After laser shut-off, generation of reverse fields has 
been observed at the front of the expanding plasma for 
ambient pressures above a certain critical value [Refs. 3 
and 21] . 

This project contains the study of field reversal in 
an expanding laser plasma. 

The spontaneous and reversed fields were detected using 
magnetic probes (inductive loops) . The probe coils con- 
sisted of five turns of #40 FormVar copper wire. 

The plasma density profiles were diagnosed with elec- 
tric double probes. The double probes were maintained at 
a fixed bias of sufficient magnitude (20 volts) to insure 
collection of saturation ion current. Since the expansion 
velocity of the laser plasma was greater than the ion 
thermal velocity, the probe current was proportional to 
the product of the ion density and the plasma expansion 
velocity. The expansion velocity was obtainable from the 
probe signals so that the relative density of the plasma 
could be calculated. 

This thesis is divided into 7 more sections. Section II 
contains the previous experimental work connected with the 
field reversal. Section III discusses the basic theory of 
spontaneous magnetic fields in laser produced plasmas, field 
reversal and shock waves in a plasma. Section IV contains 
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the details of the experimental methods. The experimental 
results are presented in Section V. In Section VI, the 
experimental results are compared with the theory. Section 
VII contains the summary and conclusions. Recommendations 
for further study are contained in Section VIII. 
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II. PREVIOUS RESEARCH 



During the latter part of 1970, as part of a study of 
the interaction between a laser plasma and the ambient 
background plasma. Stamper et. al. [Ref. 24] discovered 
that magnetic fields were spontaneously generated in the 
laser plasma. They explained the generation of these 
fields in terms of thermoelectric currents associated with 
large temperature and density gradients near the target. 

The spontaneous magnetic fields of the order of a kilo- 
gauss, were observed with magnetic probes. After this dis- 
covery, McKee [Ref. 18] made extensive studies of such 
fields at the Naval Postgraduate School. During McKee's 
investigation the following observations were made: 

(1) The maximum fields detected at a fixed position 
showed a systematic dependence on the ambient gas pressure. 

(2) The fields were symmetric about the target normal. 

(3) After several hundred nanoseconds, a relatively 
small magnitude field in a direction opposite to the direc- 
tion of the initial field appeared behind the initial 
fields. This reverse field appeared sooner in time for a 
background pressure of 250 m Torr than for 0.1 or 5 m Torr . 

McKee developed a model to explain spontaneous magnetic 
fields. According to McKee's model: If a background 

plasma is present, the expanding laser plasma will couple 
to it. If the density of the background plasma is high 
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enough so that the mean free path for binary collisions 
between the laser plasma ions and the background plasma 
ions is smaller than the dimension of the laser plasma, 
the coupling will be collisional in nature. At lower 
background pressures collisionless coupling will arise 
because the presence of the spontaneously generated magnetic 
field at the front of the laser plasma will lead to coupling 
via the ion-ion two stream instability in the presence of 
a magnetic field. As a result of the coupling ambient 
plasma will be swept-up at the front of the expanding 
plasma and a density shell will form. The shell contains 
a density gradient in the z-direction (see Figure 1) ; then 
field generation can occur in the density shell if a radial 
temperature gradient arises (see Section III for details) . 

After McKee, Bird developed another model (interaction 
model) to explain the pressure dependence of spontaneous 
magnetic fields. This model attributes the increase in the 
magnitude of the initial spontaneous magnetic fields with 
increasing background gas pressure, to the interaction which 
occurs between the laser plasma and the ambient background 
plasma. The interaction increases the axial density gradi- 
ents of the front of the expanding plasma and this, in turn, 
increases the rate at which spontaneous magnetic fields 
are generated in the interaction region. 

During Bird's investigation no density shell was detected. 
Also the shell model could not explain the field reversal. 
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The following observations on reverse field generation 
have been reported by Bird: 

(1) Field reversal is coincident with a strong inter- 
action between the laser plasma and the background plasma 
as is evidenced by a steepening and slowing of the front 
of the expanding plasma. 

(2) The onset of field reversal occurs at the front 
of the expanding laser plasma and it occurs earlier in 
space and time at higher pressures. 

(3) Reverse field can be produced long after laser 
shut-off if the laser plasma encounters a barrier. 

(4) Reverse field is initially produced at the front 
of the laser plasma streaming into a background of at 
pressures above 450 m Torr before any significant pile-up 
of the laser plasma has occurred. This narrow component 
of reverse field appeared to be generated in the swept-up 
hydrogen background plasma and it was followed by the 
appearance of a broad component of reverse field as the 
plasma front steepens and slows down due to the interaction 
between the laser plasma piston and the background. 

These results suggest that the generation of the reverse 
magnetic fields can be explained by formation of a shock 
wave driven by the expanding laser plasma piston. Such a 
shock wave compresses and heats the plasma. The formation 
and the effects of the shock waves in a plasma are discussed 
in the following section. 
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III. THEORY 



A. LASER PRODUCED PLASMAS 

The production of plasma by irradiation of a solid 
target was first treated theoretically by Basov and Krokhin 
[Ref. 1], 

The plasma is created by heating, vaporization and 
ionization of the target material. The overdense plasma 
(n > n at the surface of the target expands until the 
electron density in the front drops below the critical 
density n . Strong absorption occurs when the density is 
near the critical density n c and the plasma is heated to a 
high temperature and the pressure gradient accelerates the 
plasma. After the laser shut-off, the plasma front expands 
at a constant velocity. Thus, three distinct states can be 
observed during the course of the expansion of the laser 
plasma . 

When the expansion is into an ambient background gas 
rather than into vacuum, the last stage looks quite different. 



The critical density n for the laser plasma is the 
density at which the electron plasma frequency is equal to 
the frequency of laser radiation (n = oo 2 m e e o /e^) . At elec- 
tron densities above this critical value, the laser plasma 
reflects most of the incident radiation. The critical 
density for 1.06 pm radiation is about 10“^ cm"3. 
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At this stage, for a background gas pressure greater than 
0.2 Torr, Bobin et. al. [Ref. 2] observed a luminous sheet 
ahead of the expanding plasma. This luminous sheet lasted 
more than 150 nsec after the laser pulse peak value. Its 
velocity decreased with time. The authors attributed the 
origin of the luminous sheet to a shock wave driven by the 
plasma acting as a piston, because their observed time evolu- 
tion of the luminous front agreed both qualitatively and 
quantitatively with a simple shock wave model. 

A one-dimensional, simplified model treating the expan- 
sion of the laser plasma into vacuum has been given by 
Ready [Ref. 20]. In this treatment it is assumed that: 

(1) The plasma expands in the 2 -direction. 

(2) Incident radiation is absorbed via the inverse 
Bremsstrahlung process (free-free absorption by electrons 
in the laser plasma. 

(3) The target is initially a cold fully ionized 

plasma . 

(4) Radiation losses can be neglected. 

(5) The electrons and ions have a common temperature. 

(6) Temperature is independent of position. 



The equipartition time for transfer of absorbed energy 
from electrons to ions can be approximated by 

x £ . *25 (— ) Te- 
ei m 

where m^ is the ion mass and m^ is the proton mass. For 

and Te = 100 ev , x £ . < 1 nsec. 

ei ~ 



' 3/2 (°K) [n c (cm 3 ) Z] 1 



n = 10- 
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(7) The plasma expansion velocity is of the form 



V 



dz ,z. 
dt V 



pysically, the pressure gradients within the laser plasma 
are driving the expansion. The form of the pressure 
variation can be approximated by 



P - P s (l - §> Z < Z 



where P £ is the pressure near the surface of the target, 
and Z is the distance of the plasma front from the target 
surface . 

The temperature of the plasma can be determined by 

« 

equating the rate at which the laser beam deposits energy 
in the plasma, to the rate of change of plasma thermal 
energy plus the rate at which work W is done by the plasma. 
For one dinensional plasma of area A, expanding in the z- 
direction this gives: 



1 dW 
A dz 



1 A . 

2 dt 



(n +n . )kT - -p- 
e 'i' 3z 



dz 

dt 



( 1 ) 



From the plasma equation of state: 



p = (n i + n e )kT 



( 2 ) 



one can obtain the variations of electron and ion densities 
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(3) 



“ r 'is (1 * l> 



"es * 1 * f> 



where n- and n are the ion and electron densities at the 

1 S 6S 

target surface. The work done by the plasma in the expan- 
sion increases the plasma kinetic energy so that: 



3P # dz 
3z dt 



d r l, , , ,dZv 2, 

dt [ 2 (n i m i n e m e } ( dt J ] 



(4) 



By integrating Equations (1) and (4) from zero to Z and 
the result, 

rr 

Lj 

n = a / n s ( l - |)dz = An s z/2 , 

0 

where N is the total number of ions or electrons in the 
plasma and n is the corresponding number density at the 
target surface, Ready obtained the equations 



T 



( N . m . +N m 
r li e e 

L 6K(N.+N ) 
l e 





aV 



dt" 



[24 (N.m.+N m ) ]W 
li e e 



(5) 



( 6 ) 



These two coupled equations were solved numerically by 
McKee for the conditions of this experiment. The maximum 
temperature obtained was of the order of 100 eV and this 
maximum occurred a few nanoseconds after the arrival of the 
maximum laser intensity at the target. 
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In the above treatment, it was assumed that the 
expansion of the plasma is into a vacuum. The presence 
of a background gas (pressures above 1 m Torr) influences 
the expansion dynamics. 

B. INTERACTION BETWEEN COUNTERSTREAMING PLASMAS 

The optical radiation from a dense, hot laser plasma 
is known to be rich in the far ultraviolet and near x-ray 
regions of spectrum. Line radiation and Bremsstrahlung 
radiation from a plasma of T = 100 eV is energetically 
capable of ionizing an ambient background gas. 

Koopman [Ref. 16] has investigated the photoionization 
produced from a plane copper target by an 8-Joule, 20 nsec 
laser pulse at 6943 A° from a Ruby laser. He found that 
the ambient gas was almost completely ionized to distances 
of the order of a centimeter for a background gas pressure 

of the order of 0.2 Torr and at large radii the degree of 

-2 

photoionization fell as r . Total optical energy emitted 
between 1000-2000 A° was estimated as approximately 1±0.5 
Joule. Therefore, when a laser plasma expands into a 
photoionized ambient background gas, plasma-plasma inter- 
action must be considered on the expansion dynamics. 

According to Dean [Ref. 7] during the period of the 
laser irradiation of the target (early time) , the dynamics 
of the laser plasma expansion corresponds to the propagation 
of a radiation-driven detonation wave. After the end of 
the laser pulse (late time) the plasma expands as a blast 
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wave. This model assumes that at early times of the 
expansion phase, there is momentum transfer between the 
expanding front and the ambient background plasma and that 
as the front propagates into the background plasma, further 
energy from the laser is absorbed behind the front in which 
the density is near the critical density and is transported 
to the front by electron heat conduction. After the end 
of the laser pulse, all the momentum in the driving laser 
plasma piston is assumed transferred to the swept-up back- 
ground plasma which propagates as a blast wave. 

Hall [Ref. 12] investigated the dynamics of the expan- 
sion of the laser plasma into an ambient background of 
argon. During his investigation, a 0.5 Joule, 20 nsec Ruby 
laser was used. The target was a Tantalum slab. Using an 
image-converter camera and double probes for diagnostics, 
he found that the late time expansion dynamics of the laser 
plasma could be explained by blast wave theory. 

Koopman [Ref. 16] has suggested that the observed momen- 
tum coupling is the result of multi-coulomb collisions 
occurring between particles of the two plasmas. If a strong 
momentum transfer interaction occurs, one would expect 
"snowplowing" or displacement of the ambient plasma by the 
leading edge of the laser plasma piston, with a corres- 
ponding pile-up of piston plasma behind the interface between 
the laser plasma and the background plasma. The ambient 
and piston plasmas would be intermixed in the region of 
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interface, that is, the model is one in which the background 
plasma is "captured" at the front of the expanding plasma, 
increasing the total mass of the streaming plasma and, by 
conservation of momentum decreasing the expansion velocity. 

A one dimensional expansion model can be justified on 
the basis of the following considerations. At the target 
surface the laser beam diameter is about 1.5 mm and the 
laser radiation is absorbed within a distance of the order 
of 0.01-0.1 mm from the target [Ref. 20]. Therefore, the 
ablated target material initially will have larger dimen- 
sions in the radial direction than in the axial direction. 
Since the expansion of the laser plasma is driven by pressure 
gradients ; 



| (VP) I - -I- and I (VP) I - 

1 z 1 o ' r 1 6 

z r 

from the above discussion, 6 >> 5 , therefore 

r z 

I (VP) I >> I (VP) | where 5 and 5 are character is tic 
radial and axial scale distances. From the above inequality, 
a one dimensional laser plasma expansion can be approximated. 

When the laser plasma expands into an ambient background 
plasma, the axial position of the expanding front will be 
governed by the interaction between the two plasmas . An 
empirical expression for the late time (after laser shut- 
off) dependence of the axial position of the plasma front 
on time and on the ambient (nitrogen) plasma density, p , 
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is given by Bird [Ref. 3] : 



zcf , 0 - 1 ^- 5 

M o 



( 7 ) 



C. INITIAL PRODUCTION OF SPONTANEOUS MAGNETIC FIELDS 

In the presence of electric and magnetic fields, the 
equation of motion for the electron fluid is given by 
generalized Ohm's law. 



3V 

V"e l 3ir + <V ,,V e ! 



-en [E + V XB ] - V • P +P . (8) 

e e e ei 



where V (r,t) is the fluid velocity of electrons, E(r,t) 

and B(r,t) are the averaged electric and magnetic fields, 

P e is the kinetic stress tensor and P ^ represents the 

collisional momentum transfer between electrons and ions 

per unit volume and unit time. 

If one assumes that the electron-electron collision time 

t is much less than the characteristic time for field 
ee 

production which is of the order of duration of the laser 
pulse, then V*P e reduces to the electron pressure gradient, 
VP e - If the macroscopic flow velocity of the electron 
fluid element is small compared to the random thermal 
velocity of the electrons then the (V e *V)V e term can be 
neglected. For the aluminum laser plasma, the expansion 

7 

velocity of the front is of the order of -10 cm/sec 
[Ref. 5]. If one assumes an electron temperature of 100 eV, 
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then 




[kT /m ] 1/2 
e e 



0.01 . 



Therefore the above approximation is justified. 

If the frequencies of interest are much less than the 
electron plasma frequency, i.e., w << the electrons 

can be considered massless and react instantaneously to 
applied fields. 

P . can be approximated by the resistive drag term, 
n e eJ 01 

— — — , where a is the conductivity of plasma. Finally, 
eq. (8) can be written 



E 



— - V x b - — VP 
a e nee 

e 



(9) 



By taking the curl of E and using the Maxwell equations, 
eq. (9) becomes 



-VxE = = - ±VxJ + V X (V xb) + V x -±-VP (10) 

3t o e ' n e e 

e 

Also, from the Maxwell equations, 

VXB = U D J + 4 || (ID 

c 

The ratio of the two terms on the right hand side of 
eq. (11) is of the order 
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( 12 ) 
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y J ~ u OE ' ^ 3t ; 

o o 
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where t is the duration of the laser pulse (x 5 * 10 sec) 
o o 

for all situations of interest a > 1 mho. Therefore, 
eq. (12) becomes: 
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The above inequality shows that 



_1_ 

2 

c 



9E 

3t 



^o J 



<< 



1 



i.e., the displacement current can be neglected in compari- 
son to the conduction current. In doing so, equation (11) 
becomes 



V x b = y J (13) 

o 

using eq. (13) , the first term in eq. (10) becomes 

- I Vxj = __J^ V x(7xb) = - -^r[V(V-B - V 2 B)] 

M o ^o 

(13) 



since 
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The third term in eq . (10) can be written: 



1 , kT 

Vx— =— VP = V X -i-V(n kT ) = V x (— VT + — -V n ) 

nee rises eeriee 

e e e 
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(15) 
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The density gradient is a vector; = A but the tempera- 

kT ^e 

ture is a scalar quantity; H 0 then, using the vector 



identity; 



Vx(0A) = 0VxA + (V0)xA 



(16) 



Eq. (15) becomes: 



kT 

Vx (vHI Vl lJ 

n e e e 
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putting all terms in their final form together, eq. (10) 
becomes 



H = 1T i 5 -V 2 E + vx (V e *B) + [(|vT e )x(-r 7ne)) (18) 

o e 

The first two terms on the right hand side of eq. (18) 
are the diffusion and convection terms. There is no external 
magnetic field applied, i.e., initially B = 0. Therefore 
the generation of an initial magnetic field requires that 
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the last term, the thermal source term (S) be non-zero 
due to non-parallel temperature and density gradients. 

The initial configuration of magnetic field will decay 
away in a diffusion time t [Ref. 14], 



T 



4ttctI/ 



(19) 



where L is a length, characteristic of the spatial varia- 
tion of B. For times short compared to the diffusion time 
or, in other words, when the conductivity is so large that 
the first term in eq. (18) can be neglecged, equation (18) 
can be reduced to: 



3B 
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Vn 



= V x (V xB) + (- 
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This equation expresses the fact that the rate at which the 
fields vary in a region, depends on the rate at which the 
fields are created in the region and the rate at which the 

field flows in or out of the region. 

• • Y 

For adiabatic expansion of the electrons, PV = const. 

M 

0 y 

or n kT ( ) = const, (where M is the total electron 

e e m n e 

e e , 

, 1 — y 

mass in a given volume V) or kT^rig = const, and the 

source term vanishes (VT xVn aVri xVn = 0). Thus, non- 

e e e e 

adiabatic conditions are necessary for the generation of 
magnetic fields by the thermal source term. 

If the laser energy is supplied to the electrons at a 
sufficiently high rate, then the magnetic source term is 
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non-zero in the resulting non-adiabatic state. Stamper 
[Ref. 26] analyzed the problem by means of the electron 
heat equation. 



dW. 



dt 
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where W„ is the energy supplied to the electrons per unit 
hi 

time. The first term on the right hand side of equation (21) 
dW T 



nk 

Y-l 



dT 

dt 



dt 



where is the internal energy and the 



second term on the right hand side of equation (21) 
dri dW W 

kT-rr = — tt— where W rT is the work done by the electrons, 
dt dt W 

According to Stamper, the ratio W /W_ is the measure of 

H I 

non-adiabatic conditions. Stamper approximated W H by th 
product of the absorption coefficient k. ; (cm - ''-) and the 
radiation intensity I (Watt/cm ) , while Wj was approximated 



by 



rukT 

where x^ is an average electron heating time. If 

one assumes k v is due to inverse Bremsstrahlung absorption. 



in the underdense region, that gives: 
W t £ . . — 

H _ ei i/cn 



w. 



ei 



n kT 
e 



( 22 ) 



4 — 

where T . is the electron-ion collision time, r\ is the 
ei 



ei 



2. 7 x 10 

, - 3 . 

n i (cm ) 



-2 



Te 3//2 (°K) for n- 

i 



, n 21 -3 

10 cm 



and T 



ei 



100 eV. From footnote (2) , x • 



8.1 x io 



-10 



ei 



, , a -14 

3xio sec . ^ 

^ m . 

sec., and the ratio x £ ./x . ~ 9.25 x 10 “ >> 1 , 

' ei ei ~ m 

P 
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index of refraction of plasma (n - 1) . Actually near the 

cut-off, the index of refraction goes to zero, thus the 

• • 

ratio W„/W T >> 1 and the plasma appears non-adiabatic . 

During the time the plasma is absorbing the laser radiation 

for the conditions of this investigation, the ratio W /W 

n X 

is of the order of one and the plasma appears non-adiabatic 

• • • 

But if W H << W-j. and then the plasma appears adiabatic. 

The direction of magnetic field is determined by the 
geometry of the laser plasma. At early times, the laser 
plasma, as justified before, will expand in the axial direc 
tion and create an axial density gradient. The temperature 
gradient is due to the finite radial extent of the laser 
beam and is in the negative radial direction. The cross 
product of these two gradients will generate azimuthal and 
axisymmetric magnetic fields [Fig. 6] as observed. 

D. FIELD REVERSAL 

At later times, after cessation of the laser pulse, 
the gradients become smaller and field production becomes 
insignificant. However if the laser plasma expands into 
a background plasma of sufficient density, shock heating 
at the front will produce an axial temperature gradient 
[Ref. 21] . Since a radial density gradient always exists 
in the plume-shaped laser plasma, magnetic field in a 
direction opposite to the initial field will be generated. 
The equation describing this field is: 
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(23) 



Schwirzke [Ref. 21] and Bird [Ref. 3] observed this magnetic 
field reversal several hundred nanoseconds after laser shut- 
off. The production of reverse field at the front of the 
laser plasma requires that the ratio: 




be greater than one. 

The gradients in this ratio can be approximated by: 
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where the 6's are the characteristic gradient scale lengths. 

The ratio becomes 6 *6 _/6 *6 . For a one dimensional 

zn rT zT rn 

adiabatic compression or shock compression in the z-direction, 

5^=6 and 6 _ << 6 m . The large electron heat conduc- 
zT zn zT rT ^ 

tion provides that 6 _ > 6 and thus the ratio becomes 
c rT rn 

greater than one ie. , reverse field production is possible 
under these conditions. 

The production of reverse field appears to be a conse- 
quence of the interaction between the laser plasma and the 
ambient background plasma. The interaction results in 
slowing down and a steepening of the front of the laser 
plasma (pile-up) . 

Equation (23) predicts that the field reversal will occur 
sooner and grow to a larger magnitude as the background 
pressure increases. Bird [Ref. 5] justified this prediction 
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and also observed that 



(1) The mechanism for production of reverse fields 
is "turned-on" sooner at higher background pressures. 

(2) Field reversal is coincident with a strong inter- 
action between the laser plasma and the ambient background 
as is evidenced by a steepening and slowing of the front 
of the expanding plasma. 

The electron temperature required to produce a given 
field reversal can easily be estimated by dimensional 
analysis of eq. (23) . From that equation 

K T B q *6-6 

e ^ 8 max r z ( 24 ) 

e ~ t 

where 6 and 6 are radial and axial characteristic scale 
r z 

lengths and t is the characteristic time for reverse 
magnetic filed production. Bird estimated the required 
temperature to be of the order of 1 eV for the parameters 
of his investigation which were similar to the ones of the 
present investigation. 

E. SHOCK WAVES IN A PLASMA 

If a large concentration of energy in a continuous 
medium is suddenly released, it will spread into the medium 
and at its forefront will be a shock wave. A shock wave is 
a sudden but continuous change in the state of medium, heats, 
compresses and accelerates the medium through which it 
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propagates. The shock wave may propagate through a 
plasma. 

The normal process of shock formation is the pressure 
dependence of wave velocity, so that a finite amplitude 
disturbance with increasing pressure steepens as the higher 
pressure portions move faster than the smaller ones in 
front. In order to obtain a steady profile wave, some 
mechanism must exist to limit the steepening. This limi- 
tation on steepening can be provided by either dissipation 
or dispersion. A steady state shock occurs when a diffu- 
sion process balances the steepening. Jump conditions are 
then obtainable from the conservation of mass, momentum 
and energy. In gas dynamics, the flows are properly des- 
cribed by the Navier-Stokes equations. Any shock wave 
present in the flow must therefore also be contained in the 
solutions of the Navier-Stokes equations, under appropriate 
initial and boundary conditions [Ref. 3] . However if the 
thickness of the shock wave becomes so small as to be 
comparable to the mean free path, the Navier-Stokes equations 
become invalid in describing the flow in the shock. One 
must then use the Boltzmann equation or similar kinetic 
equations for the description of the shock structure. 

The basic features of the shock structure in a plasma 
are related to the slow character of the energy exchange 
between ions and electrons and to the high electron mobility, 
as a result of which the electron heat conduction greatly 



34 



exceeds the ion heat conduction. Maxwellian distributions 



in the electron and ion fluids are established quite 
rapidly, in a time of the order of the time between parti- 
cle collisions [Ref. 33] . On the other hand the equili- 
bration of the temperatures of ion and electron fluids 
takes place much more slowly, because of the large differ- 
ence between the electron and ion masses. This relaxation 
process determines the shock front thickness in a plasma. 

A qualitative discussion of the consequences of the 
low rate of energy exchange between electrons and ions is 
given by Zel'dovich and Raizer [Ref. 33]. In this simpli- 
fied discussion it is assumed that: 

(1) electron heat conduction does not differ from ion 
heat conduction 

(2) ionization does not take place in the shock .wave; 
the wave is propagated through a gas which is already 
ionized 

In a reference frame moving with the wave a considerable 
part of the kinetic energy of the gas entering the com- 
pression shock, is irreversibly converted into heat through 
the action of ion viscous forces . The increase in ion 

temperature across the compression shock is of the order 
m.V z 

AT^ ~ — ^ — where V is the velocity of the incident gas, 
equal to the shock front velocity and m^ is the mass of an 
ion. The thickness of the viscous shock is determined by 
the time between ion collisions, x^; it is of the order of 
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an ion-ion collision mean free path x. where v.. 

1 th 1 th 

is the thermal velocity of the ions in the compression 

shock. During the compression time x^ the ions do not 

succeed in transfering any appreciable thermal energy to 

the electrons, since the characteristic exchange time 

x . ~ (— — ) / x. is very large compared to the ion collision 

e 

time x . . 

i 

The increase in electron temperature across the shock 
due to the conversion of the kinetic energy of the elec- 
trons into heat through the action of electron viscous 

forces is negligibly small. It is of the order of 
2 

e — -r — , that is smaller by a factor of m /m. than AT.. 

K S X X 

The electron gas in the compression shock is heated by some 
other process. Since the dissipation of energy due to 
electron viscous force is very small, Zel'dovich and Raizer 
[Ref. 33] suggested that the compression and heating of 
the electron gas in the compression shock is adiabatic: 
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where T ^ and T £ are the electron temperatures in State I 

and State II, and are the electron densities in State 

I and State II and y is the ratio of specific heats. For 

p 2 

a strong shock wave, the ratio — is of the order of 4. Thus 

T p 1 

e2 x 

m ~ 2.5 for y = 1.4. Therefore in a strong shock wave 

el 

propagating through a plasma with equal ion and electron 



36 



temperatures, behind the shock there is considerable 

difference in the temperatures of the two gases. After a 

small portion has undergone the shock compression it begins 

to transfer thermal energy from the ions to the electrons 

and this leads to an equilibration of the temperatures in 

. £ 

a time of the order of the energy exchange time t The 

final temperature (= = T^) is determined by the general 

equations of conservation for shock fronts. The thickness 
of the relaxation region behind the shock in which plasma 
approaches equilibrium (T = T^ = T^) is of the order of 



Ax 



v. x e . 
1 ei 



where 




Jaffrin and Probstein [Ref. 15] investigated the structure 
of a one dimensional shock wave in a fully ionized plasma 
consisting of electrons and singly ionized ions in the 
absence of any applied external electric or magnetic fields. 

The analysis of the shock structure is based on a two- 
fluid continuum approach in which they took as a "model" 
the Navier-Stokes equation, including the electric field and 
transfer terms. In so doing the overall width of the shock 
transition region turns out to be considerably larger than 
the characteristic mean free path. They found that 

(1) for mach number M < 1.12 (weak shock) the shock 
profile is monotonic, with the shock thickness of the order 
of Z^/z w here Z ^ is the down stream (behind shock) mean free 
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m e 1/2 

path of ion-ion collisions and e = (— ) . The important 

i 

dissipative mechanisms are electron thermal conduction and 
the energy exchange between ions and electrons. 

(2) M > 1.12 {moderate to strong shocks) . There is 
still a thermal shock layer whose approximate thickness 

is l 2 / £ * In this layer the electron temperature increases 
and approaches its downstream value while the ion tempera- 
ture and velocity changes are relatively small. However, 
imbedded in this first layer there exists a second inner 
shock whose thickness is of the order of the downstream 
ion-ion collision mean free path • In this second layer, 
the dissipations are due to the ion viscosity and thermal 
conduction. The heating and compression of ions occur 
mainly in this second shock layer. However, because of 
the ions' much heavier mass, the ions do not have time to 
exchange energy with the electrons in a few mean free paths 
and the electron temperature is essentially unaffected. 

The ion and electron temperatures then equilibrate in the 
downstream continuation of the thermal layer. The thick- 
ness of this downstream layer is of the order of 

(3) Debye length of the order of the mean free path 
(X ~ l ^) . In this case spatial oscillations of wave- 
lengths which are of the order of the Debye length are 
present. The charge separation is no longer small and the 
electrical effects interact with the flow in the two 
electric shock layers. 
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The structure of a collisional weak shock in 

the presence of magnetic field which is in the plane of 

shock was studied by Hu [Ref. 13] . He found that 

the shock thickness depends on four parameters, the 

electron mean free path X , the electron larmor radius r 

e e 

based on thermal speed, the electron penetration length 
m c^ 

0 2 
r = = 5 — i-vov and the electron-ion mass ratio e . Hu 

p 2.1/2 

^ (47rn e e ) 

discovered two important features of the shock structure. 

(1) It is possible to have collisional shocks with 
oscillatory structure due to electron inertia. 

(2) It is possible to have collisionless shock with 
collisions whose thickness is smaller than the collision 
mean free path due to the complicated interactions between 
electromagnetic fields and the plasma. Often in interpreting 
experiments, shocks exhibiting these two features are con- 
sidered as collision-free shocks. 

In general, the mean free path for collisions in a 
plasma determines the type of the shock. A shock layer, 
which depends on the interaction length, can be thinner 
than a mean free path. These shocks provide examples for 
studying the mechanism of collisionless dissipation. A 
similar definition of the collisionless property was given 
by Friedman [Ref. 10]. For a given shock structure, the 
dissipation provided by classical collisions (ordinary 
momentum transfer collisions) must be less than the dissi- 
pation which is required to form a steady shock. 

The first conclusive evidence of the collisionless 
shocks come from astrophysics. The plasma flow from sun, 
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the solar wind, is highly supersonic when encountering the 
magnetic field of the earth, and a shock wave must be 
formed. However, the density of this plasma is very low 
so that the mean free path for collisions is much longer 
than the dimensions of the obstacle itself, the earth's 
magnetosphere. Therefore a classical collisional shock 
wave cannot exist. Satellite observations showed the 
presence of a well-defined shock transition of the order 
of 100 Km. 

In the theory of shock waves, it is important to determine 
the properties of the medium after passage of the shock 
in terms of the unperturbed quantities before the shock 
such as densities, pressure and temperature. Assuming that 
the shock profile is (quasi-) stationary in a certain 
reference frame, a number of relations between upstream 
and downstream quantities are given by jump conditions, 
which derive directly from the conservation of mass, momen- 
tum and energy. In usual gas dynamics and MHD-theory, 
these jump conditions are sufficient to determine the basic 
properties of the state behind the shock although a theory 
of the shock structure and the dissipation process leading 
to entropy production requires a kinetic description. 

In case of a collision-free shock wave the plasma is, 
in general, not in local thermal equilibrium behind the 
shock transition. Electron and ion pressures may be quite 
different from each other, depending on the dissipation 
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mechanism. Here, the jump conditions do not completely 
describe the state behind the shock. A dissipation of the 
shock structure itself being needed. This would, e.g. 
provide information on the ion or electron thermalization 
length. 

Plasma momentum transfer calculated to be collision 
free was observed only when a magnetic field was present 
by Friedman and Dawson [Ref. 11] during plasma gun studies. 

Dean et. al. [Ref. 8] performed an experiment which was 
similar to this present one. A Nd-doped glass laser 
(8J, 30 nsec) was used. The target was 0.25 mm Lucite 
(Cg^C^) fiber suspended in an ambient nitrogen gas. They 
found coliisionless interactions even in the absence of 
an applied field and explained this apparent discrepancy 
by the fact that spontaneously generated fields always 
existed in their experiment. They measured the thickness 
of the interaction front as approximately 1 mm and esti- 
mated the binary-collision momentum transfer mean free 
path, from a screened coulomb potential of the order of a 
centimeter. 

Shock thicknesses in collisional plasmas are usually 
taken as one-to-several collisional mean free paths. The 
difficulty of establishing a given experimental situation 
as collision free has been emphasized by Wright [Ref. 32] 
who suggested that the results of [Ref. 8] can be explained 
by collisional coupling alone. Dean et. al. [Ref. 9], 
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after Wright's suggestion reperformed their experiments. 

They changed the target to pure carbon filaments to find 
out whether the multicomponent nature of the former target 
(C^HgO^) was playing a role. As a result they found the H 
and 0 in the target was inconsequential. They also changed 
background gas and used Helium instead of Nitrogen. The 
basic phenomena were qualitatively unchanged. They sug- 
gested that the momentum coupling region which was observed, 
was not a shock. It was merely a region of momentum 
coupling between interstreaming plasmas. The mean free 
path is the e-folding distance for momentum coupling, so that 
several mean free paths are required for a high degree of coupling. 

The expanding plasma heats and transfers momentum 
to the photoionized background plasma. These inter- 
actions can be either collisional or collision-free depending 
on densities and velocities. During an experiment per- 
formed by Stamper et. al. [Ref. 27], the background densi- 
ties have typically been such that electron processes are 
collisional while the ion processes are collision-free. They 
observed a sharply defined momentum coupling interaction 
region between the laser plasma and background plasma, but 
they have not observed the existence of a shock. 

Koopman and Tidman [Ref. 17] observed a thin collision- 
less electrostatic shock front propagating radially at the 
leading edge of the plasma ball. In their investigation, 
granules of LiD have been irradiated by 15-nsec, 3-J pulses 
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from a Ruby laser. The target consisted of a cluster of 
5-to-25 ym diameter solid particles which were injected 
into the 400-ym diameter laser focal region by an electro- 
mechanical apparatus. No magnetic fields were present. 

The phenomenon is analogous to a blast wave; the expanding 
plasma acts like a piston and sweeps up the background 
plasma. This swept-up plasma gives rise to a spherical 
ion-wave shock. 

The distinction between collision-free and collisional 
mechanisms in counterstreaming plasma studies, and the 
appropriate techniques for calculating collisional cross 
sections for momentum transfer, have been the subject of 
considerable controversy. The question of the role of 
binary collisions is always a very complicated one in plasma 
experiments. There is still a need for more sophisticated 
treatments of collisional effects in plasmas. 
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IV. EXPERIMENTAL METHODS 



A. LASER 

The laser used in this investigation was a Korad-K-1500 
Q-switched neodynium-doped glass laser. The laser consisted 
of an oscillator and a single amplifier stage. The laser 
was Q-switched with a Pockels cell. The laser system was 
operated at a power level of 230 MW, which was obtained 
with a pulse width (full width at half-maximum) of 22 nano- 
seconds. The output of the system varied within about 10 
percent of 230 MW and no data were accepted if obtained at 
power levels outside this range. Other details of the laser 
are given by Davis [Ref. 6]. The laser focal spot size is 
approximately 0.15 cm in diameter after focusing with a 
30 cm focal length lens. 

B. LASER MONITORING TECHNIQUES 

The laser beam was monitored with a Korad-K-Dl photo- 
diode. This photodiode provides a signal which is propor- 
tional to the laser power and also by integrating, gives 
another signal which is proportional to the laser pulse 
energy . 

About four percent of the beam was reflected by a beam 
splitter to a MgO diffusing block for viewing by the photo- 
diode. The scattered light then passed through a 0.1 per- 
cent neutral density filter before entering the photodiode 
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as shown in Figure 2. The power signal was used to trigger 
a Tektronix 555 Dual-beam oscilloscope and was also dis- 
played on that oscilloscope. 

The energy signal was displayed on a Tektronix 564B 
storage oscilloscope. 

The photodiode energy signal was calibrated by L.L. 

McKee [Ref. 18] by using a Westinghouse RN-1 laser radiometer. 
The radiometer provided an absolute measure of the energy 
incident at the target which allowed calibration of the 
diode's integrated signal. 

The energy output of the laser was monitored on every 
shot. For energy outputs within 10 percent of 5.0 Joule, 
the pulse width was consistently 22 ± 1 nsec. A two minutes 
cooling period between shots resulted in good reproducibility. 

C. TEMPORAL REFERENCE FRAME 

The temporal relation between the diagnostic probe 
response and the leading edge of the laser pulse was deter- 
mined by taking into account the various cable delays and 
laser beam optical path lengths. A long (10 meter) delay 
line was used from the diagnostic probe to the oscillo- 
scope in order to separate the laser power signal and early 
diagnostic responses which would otherwise be superimposed 
on each other . 

By assuming that the signal travels at 1.98x10^ cm/sec 
in the RG-174/U coaxial cable, the delay was computed as 
34 1 nsec . 
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D . PLASMA CHAMBER 

The target was located in a chamber which was specially 
constructed by McKee to be used for the investigation of 
laser-plasma. The inside volume of the chamber was cylin- 
drical with a diameter of 10 inches and a height of 6 
inches (Figure 3) . 

The beam entered the chamber after passing through a 
30 cm focal length lens and hit the target at an angle of 
30° so that one can probe the plasma along the target nor- 
mal which appears to be the axis of greatest symmetry of 
the expanding plasma. 

The target used was an aluminum flat disc, 1/4 inches 
thick, and about 2 inches in diameter. It was rotated after 
a-proximately 10 shots to prevent severe cratering of the 
target in the laser impact area. 

— 6 

The chamber was pumped to a vacuum of about 2 x 10 Torr 
of air and then background gas was allowed in by a needle 
valve to fill the chamber up to the desired pressure. The 

_ g 

vacuum of about 1.3x10 Torr was obtained by using an oil 
diffusion pump. Vacuum pressures lower than 0.1m Torr 
were measured with an ionization gauge. A thermocouple 
gauge was used to measure pressures between 0.1 m Torr and 
1000 m Torr, while, above 1 Torr pressures were monitored 
by a manometer filled with diffusion pump oil. 
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E . MAGNETIC PROBES 



1 . Probe Calibration 

The magnetic probe used in this investigation con- 
sisted of five turns of #40 Formvar copper wire. The probe 
coil was about 1 mm in diameter. The probe circuit response 
was linear, with an effective area (nA) of the inductive 
loop of 5.17 5x10 m for frequencies below about 8 MHz. 

The effective area (r|A) of the probe was measured by using 
Helmholtz coils. The instruments employed for this measure- 
ment were a General Radio 1330A continuous wave oscillator, 
Hewlett-Packard 411A r.f. VTVM, Helmholtz coils and a 
Tektronix 545B oscilloscope. A circuit diagram for probe 
calibration is shown in Figure 4. 

The details of construction of Helmholtz coils and 
the use of these coils for calibration of a probe is des- 
cribed by McLaughlin [Ref. 19]. The effective probe area, 
according to McLaughlin is given by the equation 



(HA) 




( 1 ) 



where K is a parameter which differs only from coil to 

coil (for the coil which was used in this investigation, 

2 

K = 0.025 m ), E is the voltage across the probe connection 
cable and V L is the voltage applied to the Helmholtz coil. 

2 . Probe Signal Reliability and Sensitivity 

To insure that the magnetic probe signals (3B/3t) 
were meaningful, several steps were taken. The probe 
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signals were monitored under experimental conditions for 
which no target was in place, to verify that a null signal 
was obtained. Also, the probe was rotated by 180° about 
its axis and the probe signals were checked to make sure 
that they reverted in polarity. This demonstrated that 
the signals were not electrostatic in nature. Magnetic 
field signals were recorded on a Tektronix 555 dual-beam 
oscilloscope . 

3 . Probe Signal Integration 

According to electromagnetic theory, the output of 
a coil is 



e (t) 




(26) 



where r|A is the effective coil area, B is the magnetic 
field, and e(t) is the voltage generated across the coil 
which was recorded on a Tektronix 555 dual-beam oscillo- 
scope. Integration of equation (26) gives: 



B 



B 

/ dB 

0 



_ 1 _ 

nA 



t 

/ e(t)dt 

0 



(27) 



t 

The integration of the probe signal ( / e(t)dt) was 

0 

performed by measuring the area under the e(t) curve. 
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F. ELECTRIC DOUBLE PROBES 



Double probe thoery can be found in several references 
[Ref. 3] and [Ref. 4]. 

At sufficiently high bias voltage, the probe draws 
the ion saturation current I 



I 



Z e 




(28) 



where is the probe area perpendicular to the direction 
of plasma flow, n^Ze is the average ionic charge per unit 
volume in the laser plasma and is the magnitude of the 

plasma flow velocity. 

In equation (28) the assumptions are: 

(1) The areas of probe tips are equal. 

« 

(2) The plasma parameters are uniform across the probe 
tips . 

(3) Plasma-probe interactions such as secondary elec- 
tron emission can be neglected. 

The secondary electron emission coefficient (e) due to 
ion bombardment (y^) , defined as the average number of 
electrons ejected per incident ion, has been measured by 
Koopman [Ref. 16] for similar experimental conditions and 
found to fall in the range of 0.0 - 0.5. Therefore secon- 
dary electron emission will contribute to the total current 
drawn between the probe tips. After taking secondary 
electron emission into account equation (28) becomes 
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(29) 



I - (1+e) A x ni Ze V Lp 



The probe was biased at 20 volts to collect a saturated 
ion current. A large secondary emission coefficient of 
e = 0.5 was used in obtaining a quantitative idea of plasma 
density . 

The probe circuit is shown in Figure (5). The probe 
circuitry was enclosed by an aluminum housing attached to 
the probe and the entire system was isolated from ground. 
The diameter of the cylindrical tips was about 1 mil, they 
had a length of 1 mm and were about 1 mm apart. 

The double probe signals were recorded on a Tektronix 
555 dual-beam oscilloscope. To obtain an estimate of the 
plasma density, equation (t) is solved for and evaluated 
with Z = 1.5, £ = 0.5, A ± = 2.54 xlO -4 cm 2 , I = | , R = 1 fl 



= 1.09 x 10 



22 



x 



V 



V 



LP 



(volts ) 
(cm/sec) 



(30) 



where V is the voltage across a 1 ohm-resistor (see Figure 5) . 

G. DATA COLLECTION AND ERROR ESTIMATES 

The magnetic probe was inserted into the vacuum chamber 
through a small pipe which was soldered onto a brass base 
plate. 

The base plate could be moved about on an "0" ring seal 
located on the top of the chamber without breaking the 
vacuum seal. The base plate was scribed with reference lines. 
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These lines allowed positioning of the plate over a one 
millimeter grid attached to the top of the chamber. The 
position of the probe relative to the target was measured 
internally at the beginning of each data session to insure 
that the probe had not moved in its holder. This method 
allowed positioning of the probe to within ± ~ mm of the 
stated position in relation to the target. The electric 
probe was inserted into the chamber from a side window. 

The relative measure of plasma density from point to 
point is considered good due to the high reproducibility 
of the laser pulses. It is estimated that these relative 
values are accurate to within ± 10 percent. 

The laser gave very reproducible pulses; the output 
energy was allowed to vary as much as ± 10 percent. 

The area of probe tips is considered accurate to within 
10 percent. The plasma expansion velocity is considered 
accurate to within 10 percent. 
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V. EXPERIMENTAL RESULTS 



During this investigation, hydrogen gas was used as 
the background gas. An electric probe was introduced from 
the port No. 6. (Figure 3) and the magnetic probe from the 
top window. 

To decelerate the plasma front rapidly, a glass plate 
is placed into the vacuum chamber (Figure 7) . The main 
influence of the glass plate on spontaneously produced 
magnetic field is displayed in Figure 28. This figure 
indicates that the maximum value of reverse field when the 
glass plate is in place (at z = 1.15 cm) , measured with 
the probe at r= 0.4 cm, 0 = 0° and z = 1.0 cm, is greater 
than the maximum value of reverse field when the glass 
plate is not in place (at z = 7.5 cm) . 

Figures 8 and 9 display the relation of 3B/3t to the 

— 6 

plasma density for a vacuum pressure of 2 x 10 Torr, as 
a function of time. When the glass plate is in place, the 
magnitude of 3B/3t is greater than the magnitude of 3B/3t 
when the glass plate is not in place. 

Figure 16 demonstrates that the slope of the plasma 
density curve when the glass plate is at z = 1.15 cm is not 
very different from the slope of the plasma density curve 
when the glass plate is at z = 7 . 5 cm . 

Figure 17 shows that when the background gas pressure 
is 0.5 Torr (J^) / the existence of the glass plate at 
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z = 1.15 cm causes a noticeable steepening and a sudden 
increase in plasma density. 

Figures 10 and 11 show the relation of 3B/3t to the 
plasma density at a background pressure of 0.5 Torr (I^) 
as a function of time for two locations of the glass plate. 

When the background gas pressure is 5.0 Torr (E^) , the 
sudden increase in plasma density when the glass plate is 
in place (Figure 18) is much more pronounced (see also 
Figures 13, 15 and 28). For 5.0 Torr (t^) background 
pressure, the difference between the maximum value of 
reverse field when the glass plate is at z = 1.15 cm and the 
maximum value of reverse field when the glass plate is 
at z = 7.5 cm is greater than the difference between the 
same quantities for a background pressure of 0.5 Torr (I^) • 

Figures 12 and 14 display the relation of 3B/3t to the 
plasma density for 5.0 Torr background gas pressure at the 
point r = 0.4 cm, 6=0° and z = 1.0 cm for two different 
locational of the glass plate (note that each curve in this 
figure was multiplied by a different factor) . 

Figures 19, 21 and 23 demonstrate the comparison of two 
plasma density curves each for different background gas 
pressure, at the point r = 0.4 cm, 9=0° and z = 1.0 cm 
with the glass plate at z = 1.15 cm. The slope of each of 
these two curves is almost the same. The comparisons of 
the same curves when the glass plate is at z = 7.5 cm are 
displayed in Figures 20, 22 and 24. From these figures, one 
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can notice that the slopes of the two density curves each 
one for different background gas pressures are different 
from each other. 

Figures 25, 26 and 27 illustrate the comparison of 

the magnetic probe signals (3B/3t) when the glass plate 

is located at z = 1.15 cm with the magnetic probe signals 

(3B/3t) when the glass plate is located at z = 7.5 cm for 
— 6 

2x10 Torr, 0.5 Torr (E^) and 5.0 Torr (E^) background 
pressures. From these figures, one can see that the reverse 
field maximum is the greatest when the background gas 
pressure is 5.0 Torr (H^) and glass plate is at z = 1.15 cm. 

Figures 29, 31 and 33 illustrate the comparison of the 
two magnetic probe signals (3B/3t), each one for a differ- 
ent background pressure, as a function of time at the point 

< 

r = 0.4 cm, 0=0° and z = 1.0 cm and the glass plate is 
in place. An investigation of these figures indicates that 
the higher the background gas pressure, the greater the 
area under the negative portion of the 3B/3t curve, while 
the area under the positive portion of 3B/3t curve becomes 
smaller. This means that as the background gas pressure 
is increased, the reverse field maximum becomes greater when 
the glass plate is in place. Figures 30, 32 and 34 illus- 
trate the same comparison when the glass plate is at 
z - 7.5 cm . 

__ g 

For a vacuum pressure of 2 x 10 Torr, the magnitude 
of the magnetic probe signal (3B/3t) when the glass plate 
is in place is about 4 times greater than the magnitude of 
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the magnetic probe signal (3B/3t) when the glass plate is 
not in place (see Figure 25) . (Note that each curve in 
this figure was multiplied by a different factor.) 

When the glass plate is in place, the magnitude of the 

_ g 

magnetic probe signal (3B/3t) at a vacuum pressure of 2x10 
Torr is greater than the magnitude of the magnetic probe 
signals (3B/3t) at 0.5 Torr (I^) background pressure and 
at 5.0 Torr (I^) background pressure (see Figures 29 and 
31) . But when the glass plate is not in place, the magnitude 

_ g 

of the magnetic probe signal (3B/3t) at 2x10 Torr vacuum 
pressure is smaller than the magnitude of the magnetic 
probe signals (3B/3t) at 0.5 Torr (F^) background pressure 
and at 5.0 Torr (E^) background pressure (see Figures 30 
and 32). From Figure 34, it can be seen that when the glass 
plate is not in place, the magnitude of the magnetic probe 
signal (3B/3t) at 0.5 Torr (H^) background pressure is 
greater than the magnitude of the magnetic probe signal 
(3B/3t) at 5.0 Torr (E^) background pressure. (Note that 
each curve in this figure was multiplied by a different 
factor.) But Figure 33 illustrates that when the glass 
plate is in place the magnitude of the magnetic probe signal 
(3B/3t) at a background pressure of 0.5 Torr (I^) , between 
0 and 100 nsec is greater but between 100 and 450 nsec it 
is smaller than the magnitude of the magnetic probe signal 
(3B/3t) at a background pressure of 5.0 Torr (I^) • (Note 
that each curve in this figure was multiplied by a different 
factor . ) 
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Figure 35 illustrates that when the glass plate is 
in place, for a background pressure of 5.0 Torr (J^)/ the 
slope of the plasma density curve at the point r = 0.4 cm, 
9=0° and z = 1.0 cm is greater than the slope of the 
plasma density curve at the point r = 0.8 cm, 0=0° and 
z - 1.0 cm. But when the glass plate is not in place 
Figure 36) the slopes of two plasma density curves are the 
same . 

During this investigation, for a background pressure of 
5.0 Torr (I^)/ the maximum value of the magnitude of the 
reverse magnetic field was measured as 150 gauss when the 
glass plate is in place at the point r = 0.4 cm, 9=0° 
and z = 1.0 cm and at t = 440 nsec, while at same background 
pressure, point and time but when the glass plate is not 
in place, the maximum value of the magnitude of the reverse 
magnetic field was measured as 45 gauss only. 
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VI. COMPARISON OF THEORY AND EXPERIMENTAL RESULTS 



As mentioned in Section III-C, the production of initial 
field requires a density gradient in the - z direction and 
a temperature gradient in the -r direction (see Figure 6) . 

The cross product of these two vectors generates a magnetic 
field in the 6 direction i.e., in the clockwise direction 
when looking in the direction of the laser beam. 

If one of the gradients is reversed to the opposite 
direction, a magnetic field, in the counterclockwise 
direction is produced. The radial temperature gradient is 
always in the -r direction and its direction remains unchanged, 
either the glass plate is in place or not. Therefore the 
only possible gradient which can be reversed by placing 
the glass plate at z = 1.15 cm, is the axial density gradi- 
ent. The equation of continuity is 



fjr + V • (nV) 



0 



(31) 



if cylindrical symmetry is assumed; 



3_n 3_n . dz _3_n_ . dr 

3 1 3z ’ dt 3r dt 



0 



(32) 



Also for the laser-produced plasma expanding in the z- 



direction 



dz 

dt 



>> 



dr 

dt 



thus: ft ~ - ft * V where v ■‘■ s t ^ ie expansion velocity of 

the plasma front. 

Figure 13 demonstrates the relation of Bq to the 
plasma density as a function of time at the point r = 0.4 cm, 
0=0° and z = 1.0 cm for a background gas pressure of 
5.0 Torr (f^) when the glass plate is at z = 1.15 cm. 

Between t ~ 100 nsec and t = 180 nsec. ~ is positive, 
which means that is negative or the axial density 
gradient is in the -z direction . Figure 28 illustrates 
that at t = 150 nsec the reverse field appears when the 
glass plate is in place. Since the axial density gradient 
at this time is in the -z direction, the production of the 
reverse field can not be correlated to any change of sign 
in the direction of the axial density gradient, since 
there is no such change when the reverse field is ''turned 
on . 

There always exists a radial density gradient in the -r 
direction, no matter where the glass plate is (see Figures 
35 and 36). Therefore it requires an axial temperature 
gradient in the -z direction. Such an axial temperature 
gradient can be generated by the shock heating of the 
plasma . 

Figure 28 demonstrates that the magnitude of the 
initial magnetic field, when the glass plate is in place 
is smaller than the magnitude of the initial magnetic 
field when the glass plate is not in place. This decrease 
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in magnitude of the initial magnetic field is caused by 
the decrease in axial density gradient due to pile-up of 
laser plasma in front of the glass plate. This piled 
up plasma then diffuses radially and also causes a decrease 
in radial density gradient. This decrease in radial den- 
sity gradient can be seen if one compares Figure 35 and 
Figure 36. 

Figure 18 displays the steepening in the plasma den- 
sity between 150 and 200 nsec, when the expanding plasma 
impinges on a glass plate. This steepening can be related 
to the formation of a shock wave. Such a shock wave causes 
an increase in plasma temperature (see Appendix A) and this 
increase results in an axial temperature gradient in the 
-z direction, since the shock wave travels in the +z direction. 

As long as the ratio. 



is greater than one, the reverse field only reduces the 
magnitude of the initial magnetic field. But when this 
ratio decreases to a value which is less than one, the 
reverse field becomes the dominant field. 

From Figure 18, the ratio of two densities, one at 
t = 159.37 nsec the other at t = 180.83 nsec, is: 



l 

3 1 






1.5 



and 



At 



21.46 nsec 
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in order to obtain this density ratio, the shock wave 

must have a Mach number of the order of 1.37. According 

to the classification of Jaffrin and Probstein [Ref. 15] 

this shock is a "moderate" shock whose thickness is of 

the order of Z/e, where Z is the downstream (before shock) 

m e 1/2 

ion- ion collision mean free path, and e = ( — ) . The 

i 

electron temperature increases in this layer but the ion 
temperature increases in the second inner layer (see Section 
III-E) whose approximate thickness is the downstream ion- 
ion collision mean free path Z. 

For this Mach number M = 1.37, a temperature ratio, 



T . + T:, 
e2 12 

T , + T, . 
el ll 



= 1.368 



is obtained from the Rankine-Hugoniot jump conditions 
(see Appendix A, equation 3) . 

Since the shock thickness can be defined from the 
electric probe data as 5 = Vt, where V is the velocity of 
the expanding front and x is the rise time of the probe 
signal . 

If no interactions were occurring between the inter- 
streaming plasmas, then the plasma expansion velocity 
would be independent of background gas pressure. The 
observed dependence of expansion velocity can be seen in 
Figures 20, 22, 24, between 100 and 200 nsec. From these 
figures one can conclude that as background gas pressure 
increases, the expansion velocity decreases. 



60 



To have an idea about the structure of the shock 



one must compare the shock thickness with the collison 
mean free path length. The ion-ion collision mean free 
path can be calculated from 

10 3-3 2 

l (cm) s x t. (eV) , 

n i (cm '“) x in A 

where T^ is the ion temperature. In order to calculate l, 
one needs to know n . and T . . 

'i l 

1 2 

From kT^ = -^m^V > where V is the expansion velocity 
and m^ is the ion mass, T^ can be estimated. If one 
assumes that the Shock wave propagates in the background 
gas, or background gas flows into the shock relative to 

the frame of the shock, for hydrogen background gas, 

< 

-24 7 

m^ = m = 1.673 x 10 g and V t 10 cm/sec. For these 
values of m^ and V, T^ = 52 eV is obtained. Assuming a 
fully dissociated and fully ionized background gas at 
5.0 Torr the downstream ion density is 

£ 2 x 3 . 3 x 5 x 10 16 cm" 3 £ 3.3 x 10 17 cm" 3 . 

For these values of n and T^ , % is of the order of 
8.2 x io 3 cm. 

The shock thickness can be obtained from Figure 18, 

-9 7 

in which t = At = 21.46 x io sec and for V ~ 10 cm/sec 

-2 

6 = xV = 21.46 x io =0.21 cm. Since the shock thickness 6 
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is greater than the ion-ion collision mean free path, 

the formation of a collisional shock wave must be expected. 

In summary it can be said that whenever the background 
gas pressure is high enough, the expanding laser plasma 
acts like a piston and sweeps up the photoionized back- 
ground plasma. In front of this piston a shock is formed 
to drive an axial temperature gradient which is large 
enough to produce a reverse magnetic field. 
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VII. SUMMARY AND CONCLUSIONS 



The observed dependence of the plasma front expansion 
velocity on background gas pressure indicates the exis- 
tence of an interaction between the laser plasma and the 
photoionized background plasma. As a result of this inter- 
action, a collisional shock wave, driven by the laser 
plasma piston is formed. This shock wave heats and com- 
presses the plasma. 

Initially the spontaneously produced magnetic field is 
in the 0-direction (see Figure 1) and the thermal source 
term for these magnetic fields consists of the cross 
product of the radial density gradient and the axial den- 
sity gradient. When the shock wave is formed an axial 
temperature gradient is produced. This axial temperature 
gradient and the radial density gradient which always 
exists in the plume shaped laser plasma are the main cause 
of the reversal of these initial magnetic fields. 

The comparison of the shock thickness and the colli- 
sion mean free path determines the type of the shock. If 
the shock thickness is greater than the collision mean 
free path, the shock is collisional in nature. 

Determining the mean free path length in a plasma is 
always a complicated undertaking and the results should be 
considered as crude estimates. More sophisticated treat- 
ments are necessary in order to determine the collisional 
effects in plasmas. 
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The reversal of the spontaneous magnetic field indi- 
cates the existence of a shock wave propagating through 



the plasma, 
waves can be 



Such generation of magnetic fields by shock 
a very important mechanism in astrophysics. 
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VIII. RECOMMENDATION FOR FURTHER STUDY 



The electron temperature has been evaluated indirectly 
from the shock parameters. An experimental measurement 
of the electron temperature gradient is required for a 
further study of field reversal in shock waves. Plasma 
temperature, before and after the plasma front encounters 
the glass plate can be measured by means of x-ray diagnostic 
techniques or by other means. Also, more accurate data 
about the r dependence of the plasma density are needed 
to provide a quantitative comparison between experiment 
and theory. The data about the r dependence of the plasma 
density can be obtained by use of the electric probe 
technique . 
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APPENDIX A 



Rankine-Hugoniot jump conditions can be determined by 

using the equation of continuity, conservation of momentum 

density and conservation of energy density in the plasma 

for steady state conditions [Ref. 29] . 

c ^ 

The Mach number is defined for -2- = v = — as: 

c v 3 



M 



V. 



, 5k(T el +T il\ 
1 3M J 



( 1 ) 



where is the downstream (before shock) fluid velcoity, 
T ^ is the downstream electron temperature, is the 

downstream ion temperature and M is the ion mass. 

I 

The jump conditions^ are: 




+ 11 



( 2 ) 



T i2 

T il 



+T 



e2 



+T 



el 



( M 2 + 3 ) ( 5 M 2 - 1 ) 

16 m 



(3) 



where the subscript "2" implies upstream (after shock) 
value . 
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FIGURES 




Figure 1. Cylindrical polar coordinate system 
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A - Nd Amplifier Laser 
B - Beam Expansion Optics 
C - Nd Oscillator Laser 
D - Polarizer 
E - Pockels Cell 




Figure 2. Block diagram of experimental layout 
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X 




Figure 3. Top view of vacuum chamber. Port #1 is the laser 
beam entry port, #2 is the reflected laser beam 
observation port, #3 is the transmitted laser 
beam observation port, #4 is the optical observa- 
tion port, #5 is the glass plate entry port, 

#6 is the electric probe entry port, #7 is the 
probe observation port. 
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Figure 5. Electric double probe circuit 
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Target Material 





Figure 6. Geometry of a laser produced plasma and 
the self generated magnetic field 
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Electric Probe 





Expanding Plasma 



b) Side View 



Target 



FIGURE 7. Arrangement of the target, glass plate 
and probes 
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Relation of 3B/3t to the 
plasma density as a 
function of time at the 
point r = 0.4 cm, 0=0° 
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